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Abstract 

Spin-polarized band calculations for LaSr7B48 show a weak ferro-magnetic state. 
This is despite a low density-of-states (DOS) and a low Stoner factor. The reason 
for the magnetic state is found to be associated with a gain in potential energy in 
addition to the exchange energy, as a spin-splitting is imposed. An impurity like 
La DOS is essential for this effect. It makes a correction to the Stoner factor, and 
provides an explanation of the recently observed weak ferro-magnetism in doped 
hexaborides. 



Recent experimental works on doped hexaborides show many surprises |T|, 0. Among 
the unexpected properties is the observation of a very weak ferromagnetic state that 
persists up to large temperatures when small amounts of Ca are replaced by La in the 
hexaboride system, Laa;Ca(i_2,)B6 0. This is surprising since the density-of-states at the 
Fermi energy, N{Ef), and the Stoner factor are expected to be small in this system. 
Band calculations 0] in the local density approximation (LDA) ^ on SrBg and CaBg 
confirm that the DOS is low, with Ep falling in a valley of the DOS, and spin-polarized 
calculations for EuBg show a very weak spill-over of the spin-density from the Eu-f spin 
to the rest of the valence band. Electronic structure calculations based on the local 
(spin) density approximation, either by spin-polarized calculations or by calculations of 
the Stoner factor S, are in general able to detect if a system is ferromagnetic or not 
0, 1^, If. Because the Stoner factor for La3;Sr(i_j^^)B6 is far below the limit for magnetism, 
it is tempting to suggest that the observed magnetism is not a band magnetism, but 
something else, perhaps excitonic ^ or a polarization of the dilute electron-gas pi| . 
However, for Ce it is found that the Stoner factor underestimates the tendency towards 
magnetism, because the stabilization energy is not only coming from exchange, and the 
Stoner criterion has to be corrected for a potential energy Even if this correction is 
small in Ce, there is nothing in principle that prevents a large correction to the Stoner 
model in other systems. In this letter we show that La doped SrBg is such a system. 
According to the traditional understanding of ferromagnetism, the exchange energy is the 
source of a magnetic transition, but here we find that the band structure is such that the 
exchange splitting triggers off charge transfers that contribute to the total gain in energy 
via the Coulomb interaction. 



1 



Electron energies are calculated by the linear muffin-tin orbital (LMTO) method using 
the LDA. The structure of SrBg is simple cubic, with one Sr in the corner, and a Bg 
octahedron in the center. The internal parameter z which determines the shortest distance 
from one of the B-atoms to a cube surface, is in these calculations 0.215 of the lattice 
constant. Additional 'empty' spheres are inserted at the cube edges, between the Sr- 
atoms, in order to account for non-sphericity of the potential in the open part of the 
structure. Calculations where made for the basic unit cell, a double cell and a large 
supercell of 8 cubes (2 cube lengths in x,y and z-direction), containing totally 56 atoms 
and 24 empty spheres. The basis set for the largest unitcell includes s,p,d for the atoms, 
and s,p for empty spheres. The small cells include also f-states for the Sr-sites. The 
occupied f-part is always very small, even for the La impurity site. The lattice constant 
is fixed at 7.93 a.u. in each calculation. No lattice relaxation around the La impurity is 
considered. The number of k-points varies from 286 to 20 irreducible k-points for the unit 
cells of 10 and 80 sites, respectively. 

The bandstructure for undoped SrBg agree well with the one in ref. |0. The param- 
agnetic (non-spinpolarized) DOS from the supercell calculation for LaSr7B48 is shown in 
fig. ^ The total DOS at Ep, N{Ep), about 48 states/cell/Ry, is very concentrated on the 
single La site, with NiaiEp) ~ 19 states/cell/Ry, cf Table 1. On the remaining 55 sites, 
each Sr has on the average 1.0 and each B roughly 0.4 states/cell/Ry, respectively. By 
comparison, a rigid-band model obtained by adding one electron to the band structure 
without the La impurity (cf. fig. [l|) gives very different values: Of the total DOS, 32 
states/cell/Ry, each Sr has 1.4 and each B 0.25 states/Ry. Clearly, La behaves a localized 
impurity site with a very large local NiaiEp), much larger than on Sr. The additional 
valence electron remains located on the L be seen in Table 2. By comparing 

the charges in the undoped and doped supercell one finds that the remaining Sr7i?48- 
spheres receive no additional charge. The difference in charge between the two systems 
is never larger than 0.01 electron per site, among these sites. The derivative of the DOS 
becomes very large, about 1400 states/cell/Ry^, out of which 700 is coming from the La 
site, compared to about 30 for each Sr site when the rigid-band model is applied to the 
undoped bandstructure. This asymmetric and uneven distribution of the DOS near Ep 
is a condition for real space charge transfers as T increases and as an exchange splitting 
is imposed. The transfer goes in both cases to the La site from the rest of the system. 

Spin-polarized calculations for pure SrBg, and for a double cell SrLaBi2, do not give 
stable moments of a significant size. The supercell calculation of LaSr7B48 using 20 k- 
points finds a small moment of 0.10 fis/cell at low temperature. The ferromagnetic 
state is about 10 mRy/supercell lower in total energy than the paramagnetic solution. If 
fewer number of k-points are used one finds that the (paramagnetic) DOS and its energy 
derivative near Ep are larger, cf Table 1. The spin moment increases to 0.24 /le/cell for 
the lowest number of k-points, with a total energy 11 mRy lower than the paramagnetic 
state. The moment is always largest on La, with more than 50 percent of the total moment 
on and near the La impurity. The moment on B-sites far from the impurity is negligible. 
This fact and the uneven distribution of local DOS show that a rigid-band picture of 
dilute La-doped hexaborides does not apply. The bands near Ep are modified by the 
spacially non-uniform perturbation of the potential around the impurity site. 

The original Stoner model ||^ can be applied together with ab-initio paramagnetic 
band results in order to study the spin susceptibility or the onset of magnetism |^, ^. An 
applied exchange splitting, ^, leads to a transfer from minority spins to majority spins, 
and a loss in kinetic energy AK = NC,"^, where the DOS at Ep is assumed to be 
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constant. The same spin transfer leads to a gain in exchange energy Ai?^. = N'^I^^'^, 
where the exchange integral, Is, can be calculated within LDA 0. If AE^ > AK, there 
will be a total gain in energy and a ferromagnetic transition can occur. By elimination 
of factors of and one obtains the Stoner criterion, S = NIs > 1. If a partial DOS 
is not constant near Ep, there will be a possibility of charge transfers as function of ^, 



and additional contributions to the total energy. As in ref. [^, it is here possible to 
identify one sub-band, the La-impurity (mainly d) band, which has a large derivative of 
its DOS, N'^^{Ef). The total DOS at Ep is N = Npa + A^,;, where A^^, is the rest of the 
DOS. By neglect of all other derivatives one can make a model of the charge transfers as 
function of ^ and ksT. From ref |TT[, we get the charge transfers Aq^ = N^N'^^^"^ /N and 
AQ = N^N'p^TT^ (ksT)'^ / 6N , respectively, where it is assumed that N'^^ ■ ksT and A'^„ ■ ^ 
are small. 

Next we calculate the Coulomb energy associated with the charge transfer, AE{T) = 
UqAQ, from self-consistent paramagnetic calculations at different T. The charge density is 
p(T, r) = J2kj fi^p T) I ^(cp where / is the Fermi-Dirac function, and '^{e'^, r) 

is the wavefunction for point k and band j. The electronic total energy, Etot-, contains 
kinetic. Coulomb and exchange- correlation terms: 

Etot{T) = j feN{e)de - J {^V^ + e - fi)pd\ + Em (1) 

Ve is the electronic part of the Coulomb potential, Em is the Madelung part of the electro- 
static energy, and e and fi are the exchange- correlation energy and potential, respectively. 

The first term in Eq. |1], the kinetic term, becomes Ek = \{kBT)'^N{Ep), giving 
the usual Sommerfeld term in the specific heat Cei = = 7T. The difference, E{T), 
between the fully calculated Etot in eq. and Ek defines of the non-kinetic contribution 
due to charge transfers: 

E{T) = j feN{e)de - E^T) = J{^V, + e- p)pdh - Em{T) (2) 

and 

AE{T) = E{T) - E(0) (3) 

Calculations of E{T) need to be selfconsistent because of relaxation, while E^ and 
the Stoner product are only moderately affected by non-linear response effects from the 
system. At low T is AE nearly proportional to AQ, AE ^ UqAQ, while it saturates for 
higher T. The calculated parameter Uq varies between 1 and 2.5 Ry/el./La for fc^T in 
the interval 0.5 to 15 mRy. This is from calculations using 10 irreducible k-points. It is 
probable that the main part of AE is coming from the Madelung term, because of the 
sizeable charge transfer to La. This value of Uq is calculated witout taking into account 
lattice expansion, but the mechanism of large charge transfers is expected to have an 
influence on thermal expansion and specific heat. 

The charge transfer towards La is a result of the increasing La-DOS near Ep (fig. |l]). 
A very similar charge transfer pattern will result from a spin splitting, ^, of the para- 
magnetic bands. The similarity between the two types of charge transfer permits us to 
write AE^ = UqAq^ = UqN'i^^N^^"^ /N . With this energy included, one obtains a modified 
Stoner criterion ||11|| : 

Su = Nh + UqN'l,N,/N^ > 1 (4) 
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The calculated Stoner product Nig is 0.21, A?"^^ ~ 730 per Ry^- cell, N ~ 48 and 
~ 29 (per Ry ■ cell), to give Aq^ 440^^ ^nd AQ = 720(A;bT)2, with ^ and ksT in 
Ry. The values make the Coulomb term in eq. § very large, about 14. However, even 
if the Coulomb energy represents a major correction in this case, there are several facts 
indicating that other energies are involved as well. First, the exchange splitting ^ is not 
uniform over all sites. The spin-polarized calculation shows that ^ is finite only at or very 
near to La, and the energy of a variable ^ is not counted. Second, although the amount 
of charge transfers due to ^ and ksT are similar, there are long tails of the Fermi-Dirac 
distribution which make the radial dependencies of Aq^ and AQ{T) different. Third, the 
effective derivative of the La-DOS, N'i^^{e,T), (taken as /(— ^)A^j(e)(ie), is likely to be 
largest at low T, because the bare DOS has not the same large derivative over a wide 
energy range. The mechanism of charge transfer and magnetism is therefore expected 
to disappear at some large T. This is confirmed in the spinpolarized results. By raising 
ksT to 3 mRy in the calculation with 10 k-points, the magnetic moment goes to zero. 
It is also interesting to note that the value of N'^^ is larger in the calculation with fewer 
number of k-points, as shown in Table 1. The number of k-points are insufficient for a 
good convergence in this case, but the variations of N'j^^ and the magnetic moment in the 
spin-polarized results as a function of the number of k-points, are similar. 

The dependence of the charge transfer is a result of assuming at most a linear DOS 
variation in the simple model. In reality, there are more complex variations of the DOS 
and its derivatives around Ep- By imposing a ^ on the calculated bands rigidly, one finds 
that Aq^ is proportional to only below ~ 1 mRy, it is linear around the interval 2-4 
mRy, and it saturates at about 0.007 el. /La for ^ larger than about 5 mRy. By comparing 
the energies for large ^ one can estimate the maximum for a stable moment: 

N^hi'^ + f/o0.007 = Ni'^ (5) 

With ?7o=1.5, this gives ^=16 mRy and a moment of ~ 0.8 per cell. These values 
are much larger than the values from the stable spin-polarized solution, about 3 mRy 
for ^ and 0.1 /i^ for the moment, located on La only. The model can reproduce the 
essential dependencies of the magnetic instability, and it helps to understand the reasons 
behind the transition. But for quantitative values one has to rely on the spin-polarized 
calculations. 

The magnetic state in the supercell is fragile, and in a smaller cell not even stable. 
The fact that an insufficient number of k-points can make the moment to increase, is 
not meant to represent the real physical situation. It serves to demonstrate that the 
moment depends on the derivative of the DOS, and thus to confirm the hypothesis that 
the charge transfer contributes to the stabilization of the magnetic state. The magnetism 
is in this case sensitive to other features in the electronic structure then what is normal 
for ferromagnets. A large Stoner factor needs a large DOS, usually a peak, which often 
is sensitive to disorder and thermal smearing. But a derivative of the DOS can be large 
in a fairly wide energy interval, aside a peak in the DOS. The contribution to magnetism 
from charge transfer energies may therefore resist well to thermal effects, and thereby be 
consistent with the observation of high Curie temperatures. 

When it comes to explain the properties of the real hexaboride system from the band 
results, one has to rely on extrapolations. The very low doping (x=0.01) for optimal 
magnetism in real La2:Ca(i_j;)Be 0], is much smaller than the doping in the largest super- 
cell (x=0.125) studied here. Calculations for realistically large supercells are not possible 
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at present. However, there are several indications that the mechanism favorable to mag- 
netism will remain or even be stronger at lower doping. First, by increasing the La content 
to a; = ^ in a double cell LaSrBi2, one finds a more rigid-band like La-DOS than for the 
dilute impurity in the supercell, and the local La-DOS resembles more the Sr-DOS. Effects 
of charge transfers will therefore be less pronounced than for a single La-impurity, and 
no magnetic state was found in this case. Second, the increasing La-La distance for less 
doping, will maintain a narrow, impurity-like DOS localized on La. The Fermi-energy 
will not cut this band in a rigid-band manner at a lower energy, where the DOS and its 
derivative is smaller, even though the number of electrons per average site is decreasing 
with lower doping. Instead, since the additional charge remains localized at the impurity, 
one expects that the position of Ep remains rather fixed relative to the La-band, with 
large local values of A^^^ and iV£„ even for lower doping. On the other hand, in the large 
volume far from the La-site one expects that charges and DOS conditions will be as for 
the undoped material, not knowing about the effects of La in part of the system. In that 
case will approach A^, where in N only the part of the DOS which is affected by the 
impurity should enter in eq. ^, and the correction to the Stoner factor will remain large. 

The appearance of ferromagnetism is therefore a natural consequence of the energy 
changes associated with an exchange splitting. These energies can be estimated within 
a mean-field theory. Favorable conditions for magnetism are non-constant DOS at Ep 
to make transfers possible, and a negative change in total energy as function of the 
transfers. The non-rigid band like electronic structure around the dilute La impurity 
shows up in the distribution of the ferromagnetic moment. According to these results, it 
is mainly La and the nearest sites which takes the major part of the total moment. B-sites 
remain essentially non-polarized. Alternative theories based on the excitonic mechanism 
0] need to take into account the non-rigid band effect of the La-doping. If a rigid-band 
mechanism did apply, one should expect a uniform distribution of the moment. Theories 
of spontaneous ferromagnetism of the dilute electron gas have also been proposed 0, ^ . 
However, the very large electron-gas parameter necessary for this (> 30) is much too 
large. It is the total (and not only the part due to doping) electron density that should be 
considered, and for this type of compound is smaller than 1.5 within the atomic sites 
and about 3 within the empty sites. 

What is shown here is that ferromagnetism occurs in the supercell of LaSr7B48. This 
is despite the fact that according to the low Stoner factor, ferromagnetism should not 
appear. By doing the calculations with different number of k-points it has been possible 
to show that the actual derivative of the local DOS at Ep, which is conditional for charge 
transfers, scales with the stability of the magnetic moment. This demonstrates that the 
mechanism for a ferromagnetic is not exclusively found in the exchange energy. The 
necessary condition of a localized, non-constant La-DOS near Ep, appears to be enforced 
as the doping decreases, making the mechanism of charge transfers to persist at lower 
doping. A computational verification of this will need a very large unit cell or a Green- 
function method for the limit of dilute doping. 
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Table 1: Paramagnetic DOS values (states/cell Ry, and states/cell Ry^), spin-polarized 
exchange splitting on La ^ (mRy) and magnetic moment on La and total (/i^), as function 
of the number of k-points. To the impurity site is counted also the contributions from the 
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Table 2: Differences in valence charge (el. per site) between doped and undoped SrBg- 
supcrccll for selected sites calculated from 20 k-points. To the impurity site is counted 
also the charge within the 6 closest empty spheres. 
Impurity-site (La or Sr) Sr near La Sr far from La B 

1.19 0.005 0.01 -0.005 



Figure 1: Paramagnetic, total and partial La-DOS near Ep for LaSr7B48 and total DOS 

in Sr8B48 as obtained from tetrahedron integration from 20 k-points. The zero energy is 
at Ep for LaSr7B48 and at the energy corresponding to 1 additional electron above Ep 
for Sr8B48. 
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